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The center of mass of a bright sohton in a Bose-Einstein condensate may reveal Anderson local- 
ization in the presence of a weak disorder potential. We analyze the effects of interactions between 
two bright solitons on the Anderson localization phenomenon. Perturbation calculus shows that 
even very weak interactions modify localization properties of the system eigenstates. For stronger 
interactions, i.e. when the solitons are close to each other, the localization is totally broken. It 
implies that in order to experimentally observe the Anderson localization effects, a single bright 
soliton has to be prepared and excitation of soliton trains must be avoided. 
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I. INTRODUCTION 

Fifty years ago Anderson discovered that the transport 
of non-interacting particles can be totally suppressed in 
the presence of a disorder [l[ . Waves fail to propagate in 
a disorder medium, due to destructive interference, and 
become exponentially localized. Anderson localization 
possesses different properties in different dimensions 
In one-dimension (ID) any amount of a diagonal disor- 
der leads to Anderson localization. Particles propagating 
with momentum fc in a disordered medium multiple scat- 
ter, undergo diffusive motion and finally localize with an 
exponentially decaying density profile 0-0|- Anderson 
localization is a single particle phenomenon. Interaction 
between particles can break the localization [Mil. Also 
long-range correlations of a disorder |12l-[la[or the pres- 
ence of non-linear terms in wave-equations |16l4l8| can be 
responsible for the transition from localized to extended 
states. 

Ultra-cold atomic gases are ideal systems to theo- 
retically and experimentally investigate quantum many- 
body phenomena. Anderson localization of such matter- 
waves was realized in a laboratory first in ID (T^ll^l and 
recently also in three-dimensions [2l|, In these ex- 
periments, in order to get rid of particle interactions, ei- 
ther Feshbach resonances were employed or a low atomic 
density limit was reached. Particle interactions, which 
are harmful to experiments with Anderson localization 
of atomic matter-waves, become vital to realize solitonic 
states in Bose-Einstein condensates (BEC). In ultra-cold 
atomic gases both dark and bright solitons were demon- 
strated experimentally [23l-[29| . An analysis of quantum 
character of the degree of freedom that describes position 
of the solitons allows for invention of interesting ideas for 
experiments. Bright soliton scattering on a potential bar- 
rier leads to a superposition of macroscopically distinct 
objects [sol - is^ . Interaction between two bright solitons 
allows for preparation of quantum entanglement of a pair 
of such macroscopic particles |33l | . Bright or dark solitons 
can Anderson localize in the presence of an external dis- 
order potential [sI - Ist} 



In the present article we concentrate on the Ander- 
son localization phenomena of a pair of bright solitons in 
the presence of an optical speckle potential. Two bright 
solitons interact between each other and the interaction 
can lead to delocalization of their centers of mass. The 
soliton interaction is unusual as compared to typical par- 
ticle interactions. Apart from the relative distance, the 
interaction potential also depends on the relative phase 
between solitons. 

The paper is organized as follows. In Sec. II we de- 
scribe shortly an improved version of the Bogoliubov ap- 
proach that allows for the description of quantum bright 
solitons. In Sec. Ill we present an effective quantum 
Hamiltonian for two interacting solitons in a weak dis- 
order potential. Analysis of the delocalization effects us- 
ing time-independent perturbation theory and by means 
of numerical integration of the Schrodinger equation is 
described in Sec. IV. We conclude in Sec. V. 



II. QUANTUM BRIGHT SOLITON IN A 
DISORDER POTENTIAL 

Let us consider a Bose-Einstein condensate in an effec- 
tively ID box of size L with attractive interactions be- 
tween atoms. In the mean-field description, an TV-body 
quantum state is a product ip{zi,t)ip{z2,t) . . .'ipiz]^,t), 
where iJjiz^t) fulfills the ID Gross-Pitaevskii equation 

(GPE) [allll 
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In Eq. ([T|), and in the entire paper, we have adopted the 
following units for energy, length and time. 
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where L>j±_ denotes a transverse harmonic confinement fre- 
quency, a is the atomic s-wave scattering length and m 
stands for the mass of the atoms. The ground state so- 
lution of the GPE is a bright soliton that describes a 
localized wave-packet which evolves without changing its 
shape, 
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where the chemical potential ^ = — iV^/8. We assume 
that the norm of the solution equals the total particle 
number, i.e. (V'o|'0o) = ^ ■ The soliton has two degrees 
of freedom, i.e. the center of mass position r and the 
phase 0. The bright soliton can be considered as the A^- 
body bound state where an ionization of a single particle 
costs energy \^\ (40| . 

In the Bogoliubov decription of the system there ap- 
pear two zero modes. The first is related to the global 
gauge symmetry, the other originates from the transla- 
tional symmetry of the A^-particle system [4l|, . Both 
these symmetries are broken by the classical solution ^ 
where a specific phase of the wave-function and a con- 
crete position of the soliton center are chosen. In order 
to recover the symmetries and describe degrees of free- 
dom related to the zero modes in a non-pert urbative way, 
we follow the Dziarmaga's approach [4J|. This approach 
allows us to obtain a simple Hamiltonian that describes 
the center of mass of the soliton and its phase in the 
presence of a weak external potential V{z) 
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where ipQ is the solution @ and 
mitian operator 

Pr = —idr, 
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-i/N. The Her- 
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is the momentum operator of the soliton center of mass. 
The other Hermitian operator, related to the soliton 
phase, is the particle number operator 



Ps = -ids=N-N. 
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We may restrict ourselves to the Hilbert space with ex- 
actly N particles, where = 0, because = 0. 

The external potential V{z) also couples the center 
of mass degree of freedom with the Bogoliubov quasi- 
particles. However, this coupling has been neglected be- 
cause if the strength of the potential is much smaller than 
the modulus of the chemical potential of the system, ex- 
citations of the quasi-particles can be omitted. Indeed, 



the lowest quasi-particle energy equals = and 
for |V^(2:)| <C the perturbation is not able to ionize a 
particle from the solitonic bound state [1^ [s^ . 

The eigenstates of the A^-particle system can be writ- 
ten as follows 
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where |0) b is the Bogoliubov quasi-particle vacuum state, 
^Ar(0) = 1/ ■\/27r ensures that we describe the system with 
precisely N particles, i.e. P^£,n = 0, and is an eigen- 
state of the center of mass Hamiltonian (|6]). The single 
particle density can be approximated, to the leading or- 
der in iV, by 

p{z) = (^t(z)V'(z)) « / \Mz - r)\' dr, (10) 



where ipi^) is the bosonic field operator. 

Equation ^ shows that the mass center of the soliton 
experiences a potential which is a convolution of the orig- 
inal external potential with the soliton profile. When the 
external potential is a disorder potential, e.g. an optical 
speckle potential [13, [l^ , the center of mass of the soliton 



becomes Anderson localized (3J| . The eigenstates of the 



Hamiltonian ([5]) have a shape with an overall exponential 
envelope 
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where rg is the mean position of the soliton and is 
the localization length. In the Born approximation and 
for the soliton momentum Pr > A/(27r) but smaller than 
the inverse of the correlation length of the original speckle 
potential, the localization length increases exponentially 
with Pr 
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Similar effects are also predicted in the dark soliton case 
(37| . The position of a dark soliton becomes Anderson 
localized in the presence of a disorder. 



III. QUANTUM DESCRIPTION OF A PAIR OF 
BRIGHT SOLITONS 

Let us consider a pair of identical bright solitons where 
one is located at ri and the other at r2- If I?'! — 
is much greater than the soliton size, a superposition of 
wave- functions ([4]) is a good approximation for the double 
soliton solution of the GPE ([T]). In the quantum descrip- 
tion, following a similar approach as in the single soliton 
case, we obtain the quantum Hamiltonian that describes 
the centers of mass and the phases of two independent 
solitons in the presence of a weak external potential 
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where (pi^s stand for the phases of the sohtons and r^'s 
for the positions of their mass centers. The operator is: 
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-A/N. The Hamiltonian is: 
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Hr, = ^ + J Viz)\M^-n)\^dz, (15) 

where Pr- = —idn- The independent sohtons can Ander- 
son locahze in the presence of a weak disorder potentiaL 
However, in order to analyze the locahzation properties 
of the sohtons we have to take into account their mutual 
interaction. One soliton experiences the presence of the 
other one because the tails of their profiles overlap. It 
results in an interaction potential which drops exponen- 
tially with an increase of the sohtons relative distance 
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Interestingly, the character of the interaction depends on 
the solitons relative phase. For 01 — (/)2 « they attract 
each other but if the relative phase is close to tt there 
are repulsive interactions. The Anderson localization of 
the non-interacting solitons takes place in the ri and r2 
space only because there is no external disorder in the (j>i 
and (j)2 space. We will see that because the interaction 
potential depends on the phases, the interaction induced 
delocalization is different as compared to a typical situa- 
tion where only relative particle distance determines an 
interaction potential [8l-[ll|. 

The total particle number is a constant quantity. In 
order to take advantage of this constant of motion let us 
express the Hamiltonian of the system in terms of the 
average and relative phases 
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with the corresponding momenta 

= P^, + P^, =Ni+N2- 2N, (19) 
P, . . ,20, 

The final Hamiltonian of two interacting bright solitons 
in the presence of a weak external potential reads 



where 



H = Ho + U, 
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Hi = Hr,~P^ j V{z)dN\Mz'ri)\'' dz, (23) 
H2 = Hr, +Pv I V{z) dN\Mz - r2)? dz, (24) 



and = —2/N. In the derivation of the Hamiltonian 
(l22t we have restricted it to the 2A^-particle Hilbcrt space 
where Pg = that is allowed because [Pg, H] = 0. Con- 
trary to P<f , the operator P^ does not commute with H 
and particles can be transfered between the solitons dur- 
ing the time evolution of the system. 

Having a quantum state ^'(ri, 7'2, (ys), obtained from the 
diagonalization of the Hamiltonian (PT|). we are able to 
calculate the single particle density which, in the leading 
order in N, reads 



(\Mz-n)\' + \Mz-r2)\^) 

X |\I'(ri,r2,V3)P dridr2dtp. 
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IV. 



IMPACT OF SOLITON INTERACTIONS ON 
ANDERSON LOCALIZATION 



In the presence of an external disorder potential V{z) 
the center of mass of a quantum bright soliton may re- 
veal Anderson localization. When two such solitons are 
prepared in the presence of a disorder, the localization 
effects can be destroyed due to mutual interaction be- 
tween them. In the following we analyze the impact of 
the interaction on the Anderson localization phenomenon 
within the time-independent perturbation approach and 
by means of numerical integration of the Schrodinger 
equation generated by the Hamiltonian (PT|) . 

We will consider a finite system of length L and fo- 
cus on a disorder potential corresponding to an opti- 
cal speckle potential [13, [3. Such a potential can be 
created by shining laser light on the so-called diffusive 
plate. In the far field the light intensity fluctuates in 
space which is experienced by atoms as an external dis- 
order potential. The potential is characterized by zero 
mean value V{z) = 0, where the overbar denotes an en- 
semble average over disorder realizations, standard de- 
r,,, .„i 1/2 

viation Vq 



v{zy 



and autocorrelation function 



V{7J)V{z' 2) = Vo^^^^lfp^ where ctq is the correlation 
length of the disorder. An example of an optical speckle 
potential is plotted in Fig. [TJ Interestingly, properties 
of a speckle potential can be easily modifled. Especially, 
one can create a disorder where Anderson localization 
operates as a band-pass filter that allows for the realiza- 
tion of a matter-wave analog of an optical random laser 

Eilia. 
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where 



Figure 1: Top panel: an example of an optical speckle po- 
tential V{z) for correlation length ao = 0.01 and strength 
Vb = 1.125. Bottom panel: disorder potential experienced by 
a center of mass of a bright soliton, i.e. after convolution with 
the soliton profile, K = / V{z)\'ipo{z - r)p, for N = 60. 



A. Time independent perturbation approach 

In this section we investigate eigenstates of the Hamil- 
tonian ((2T|) treating the interaction potential C/ as a small 
perturbation. The energies and eigenstates of the unper- 
turbed system ([/ = 0) read 



E A- E A- E 



(26) 



where ((p) is an eigenstate of the P^p operator related 
to the eigenvalue n^p. The energy En^ is an eigenvalue 

of the _P^/(2m<p) operator. Wc assume that the numbers 
of particles that form each of the two solitons is equal 
to N . Thus, the unperturbed system is characterized by 
= Q and '^o(v) = l/\/27r. The states Xuxn^ifi) and 
, {r2 ) are eigenstates of the Hamiltonians ([23]) and 
([M)) . respectively, and i?„i„^ and En^n^ are the corre- 
sponding eigenenergies with = 0. These eigenstates 
are Anderson localized in the presence of the disorder. 
The first order energy correction vanishes. 



Tiin20 V nin20 



= 0. 



(27) 



This is because U ^ cos 09. The only states 4'^? 

coupled by the interaction potential to ^'i^'^n^o correspond 
to = ± 1 . The lowest order energy correction is of the 

second order -E„/„20- In the lowest non-vanishing order 
the system eigenstates read 

+ E (29) 
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Reduced probability density for finding the first soliton 
localized at ri is 
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and the probability K2{t2) for finding the second soliton 
at r2 can be calculated analogously. 

As an example we choose parameters of an optical 
speckle potential like in Ref. [l^ where Anderson local- 
ization of a BEC has been realized. The length of the 
system, in the units L = 28 (i.e. 4mm), the trans- 
verse harmonic confinement corresponds to uj± = 27r x 70 
Hz and the correlation length of the disorder (Tq = 0.01 
(i.e. 1.32/im). We assume that each of the two solitons 
consists of = 60 ^^Rb atoms with the modified s-wave 
scattering length by means of a Feshbach resonance, i.e. 
a = —5.77 nm. The optical speckle potential is created 
by blue detuned light with an intensity which results in 
the disorder strength Vq = 1.125 that is 400 times smaller 
than the chemical potential of the system. 

Diagonalization of the Hamiltonians ([23])- ([M]) allows 
us to obtain unperturbed Anderson localized eigenstates 
Xnio(''i) and Xn2o('^2)- We have chosen two pairs of these 
eigenstates for further analysis. The first pair is related to 
eigenenergies Emo = 9.93 and En^o = 14.52 and the esti- 
mated localization lengths Imo ~ 0.020 and Z„20 ~ 0.026. 
They are relatively strongly localized on a scale compa- 
rable to the soliton size 2/A^ = 0.03. The other pair is 
related to i?„^o = 42 and En20 = 46.6, and /„jO ~ 0.088 
and Z„20 ~ 0.095, respectively. Both pairs of the eigen- 
states are shown in Fig. [2] and Fig. [3] The relative dis- 
tance between the solitons in each pair is much larger 
than the soliton size which implies that the soliton inter- 
action is extremely weak and the perturbation approach 

is justified. In the considered cases, the second order en- 
(2} 

ergy corrections £^„/„20 ''' orders of magnitude smaller 
than energy gap between E^^^^^q and the nearest neigh- 
boring level. 

Figures [IHSl also show the reduced probability densities 
Ki{ri), Eq. ([5^ . Dclocalization effects of the solitons are 
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hardly visible when we look at the densities in linear scale 
because the soliton interactions are very weak. However, 
logarithmic plots indicate that characteristic exponential 
profiles are not present. The mechanism of the break- 
down of Anderson localization in a finite system, con- 
sidered here, is the following. The interaction potential 
couples an unperturbed state 5'|°''„,q to states ^ = 

Xn[iXn'^i^i and *i»„»_i = X<-iXni'-i^-i where wave- 
functions Xn'.i{fi) and Xn'-'-iii^i) reveal greater local- 
ization lengths the higher the corresponding eigenener- 
gies En'i and are. If the localization lengths 

start to be comparable to the system size L the wave- 
functions X)i'i(^i) and Xn"-i(^i) become delocalized, ex- 
tended states. Contributions of the extended states in 
(PS)) are responsible for the breakdown of the exponen- 
tial decay that can be observed in Figs. [2][3l In the un- 
perturbed state '^nin20-' numbers of particles form- 
ing each soliton is well defined and equal to iV = 60. 
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Corrections, a„i„ii«'„,„, i anu ^ „.-„-._i, 

unperturbed state indicate that the dclocalization is ac- 
companied by transfers of a particle from one soliton to 
the other. By decreasing the distance between the un- 
perturbed, localized solitons we increase the interaction 
energy and breakdown of exponential profiles start closer 
to the localization centers. If the unperturbed solitons 
are located far away from each other it may happen that 
the corrections to the unperturbed state are too small to 
disturb exponential localization in a finite system. This 
is in contrast to a more typical situation where an in- 
teraction potential depends only on the relative distance 
between particles. In fact, suppose for a moment that U 
is only a function of |ri — r2|, then the degenerate un- 



perturbed states "^n^niO ^'^'^ ^n^niO ^''^^ couplcd by U. 
Diagonalization of U within the degenerate subspace re- 
sults on a linear combination of these states. Thus, even 
if C/ — >■ 0, each soliton is described by Ki{ri) which is a 
sum of the densities |xnio('"i)P and \Xn2o{i"i)\'^ that are 
localized at two different positions. In the case consid- 
ered in this publication such a mechanism is not present 
and the breakdown of the Anderson localization is an 
higher order effect. 
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Figure 2: (Color online) Top panels: eigenstates, X"io(''i) and 
Xn2o(''2), of the Hamiltonians (|23p and (|24p corresponding to 
the eigenvalues En^^o = 9.93 and E„20 ~ 14.52. Localiza- 
tion lengths are In^o ~ 0.020, ln20 ~ 0.026. Bottom panels: 
probability densities |Xnio(''i)|'^ and |xii2o(»'2)P (black, bot- 
tom lines) and the reduced probability densities for centers of 
mass of solitons Ki(ri) and K2{r2) (red, upper lines) obtained 
within the perturbation approach. The correlation length of 
the disorder potential is ao — 0.01 and its strength Vb — 1.125. 
Total number of particles in the system 27V — 120. 




Figure 3: (Color online) The same as in Fig. 2 but for En^o = 
42 and -B^jO = 46.6. Localization lengths of the eigenstates 
Xnio{ri) and Xn2o(»'2) are l„^o ~ 0.088 and l^^o ~ 0.095, 
respectively. 



B. Beyond perturbation approach: numerical 
integration of the Schrodinger equation 

The interaction induced dclocalization effects, de- 
scribed in the previous section, are very weak. This is due 
to the fact that by employing the perturbation approach 
we have to restrict ourselves to situations where unper- 
turbed solitons are localized far away from each other. 
In this section we will see that if the solitons are initially 
localized sufficiently close to each other the interactions 
can totally destroy the localization of solitons. 

In the present section we show results of numerical in- 
tegration of the Schrodinger equation corresponding to 



the Hamiltonian ([?!]). We assume that initially the sys- 
tem is prepared in the following state 

'^{ri,r2,ip) = Xnio{ri)Xn2o{r2)£.{(p), (33) 

where XiiiO are eigenstates of the Hamiltonians ((23)) - 
(P^ corresponding to eigenenergies En^o = 10.44 and 
En2Q = 12.88. They are Anderson localized around 
ri « 0.18 and r2 « with the localization lengths 
^TiiO ~ 0.029 and ^n^o ~ 0.032. The parameters of the 
disorder potential are the same as in the previous sec- 
tion, however, the size of the system has been reduced 
to L = 8. In numerical integration of the Schrodinger 
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Figure 4: (Color online) Reduced probability densities Ki(ri), 
Eq. p2|) . for different moments of time. In each panel, from 
the bottom to the top, the curves are related to: t = (black), 
t = 0.5 (red), t = 1 (green) and t = 1.25 (blue). The unit 
of time is to — 29s, see Eq. ([2}. Panels (a) and (b) show the 
reduced probability densities Ki(ri) and K2(''2), respectively, 
and are related to the case where the solitons attract each 
other, i.e. the initial state in Eq. (|33|l is chosen as the 
ground state of the effective Hamiltonian (|34|) . Panels (c) and 
(d) are related to the case where the solitons initially repel 
each other, i.e. the initial state £,{'fi) is the ground state of the 
effective Hamiltonian H34|) shifted so that the maximal density 
is located at = tt. The initial states x»iio('~i) and Xn2o(''2) 
in Eq. pSjl correspond to En-^o = 10.44 and En^o = 12.88 and 
their localization length are l„iO ~ 0.029 and Z„20 ~ 0.032. 
Size of the system is L = 8, the correlation length of the 
disorder potential ao = 0.01 and its strength Vb = 1.125. 
Total number of particles in the system 2N — 120. 




Figure 5: (Color online) Reduced probability density 
Eq. (|36|) . for different moments of time. In both panels black 
circles are related to t — 0, red squares to t = 0.5, green 
diamonds to t = 1 and blue stars to t = 1.25. Panel (a) 
corresponds to the results presented in Fig. |4^-b, i.e. to the 
case when the solitons attract each other. Panel (b) is related 
to the data shown in Fig. |4j;-d, i.e. to the case when the 
solitons initially repel each other. Note that the range of the 
horizontal axis in the panels is different, i.e. (— 7r,7r) in (a) 
and (0, 27r) in (b). 
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Figure 6: Evolution of the IPR, see Eq. Panels (a)-(b) 

correspond to the results presented in Fig. |4^-b, i.e. to the 
case when the solitons attract each other. Panels (c)-(d) are 
related to the data shown in Fig. [ll:-d, i.e. to the case when 
the solitons initially repel each other. 
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Figure 7: (Color online) Single particle density p{z), Eq. (|25[) . 
for different moments of time. In both panels red solid line 
is related to t — 0.5, green dashed line to t — 1, blue dot- 
ted/dashed line to t = 1.25. Panel (a) corresponds to the 
results presented in Fig. |4^-b, i.e. to the case when the soli- 
tons attract each other. Panel (b) is related to the data shown 
in Fig. St-d, i.e. to the case when the solitons initially repel 
each other. 



equation the space (ri ,r2,f) has been discretized and 
we have adopted absorbing boundary conditions at the 
ends of the axes and periodic boundary conditions in 
the If degree of freedom. The numerical integration is 
time consuming because in order to reproduce the corre- 
lation length of the disorder, the number of points in the 
ri spaces have to be very big. On the other hand, number 
of points corresponding to the (f degree of freedom can 
be chosen to be quite small. 

Time evolution of the relative phase f is much faster 
than the evolution of the centers of mass of the solitons 
because the corresponding mass parameters are very dif- 
ferent, i.e. m^p = ~2/N while the soliton masses are 
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equal N. We choose as an initial state (,{ip) the ground 
state of an effective Hamiltonian 

H^.eS = 77^ + UeE{^) + P^fi, (34) 

where 

P = J dzdr [|x„2o(OP - IXnio(OP] X 
xV{z)dN\Mz-r)\\ 
Ucs{^) = j dridr2 |Xnio(?'i)P lx«2o(?'2)px 

x[/(ri,r2,iy9) (X — cos(y9. (35) 

Such a ground state is strongly peaked around if ^ Q 
which implies that initially the solitons attract each 
other. In time-evolution of the state (p3|) we may expect 
that the reduced probability density 

i^d^) = j \^{ri,r2,(p)f' dridr2, (36) 

follows the probability density of an instantaneous 
ground state of the Hamiltonian ([34]) where Ucs and /3 
are obtained by substituting in Eq. ((35|) |Xnio(''i)P by 
time-evolving Ki(ri). 

Figure |3h.-b shows time-evolution of the reduced prob- 
ability densities Ki{ri). Interaction between solitons is 
responsible for the breakdown of the initial Anderson lo- 
calization of the solitons. The delocalization effects are 
much stronger than those analyzed in the previous sec- 
tion. Plots of the reduced probability density {(p) pre- 
sented in Fig. [5^, indicate that while K^{ip) does not pre- 
cisely follow the density of the instantaneous ground state 
of the initial degree of coherence between the soli- 
tons is practically preserved in the time-evolution. The 
density K^{ip) is concentrated around zero, i.e. the rela- 
tive phase between the solitons ip fa and consequently 
the solitons always attract each other. 

In Figs. I311S] we also show the evolution of the reduced 
probability densities in the case when the solitons ini- 
tially repel each other. That is, as an initial state ^((p) 
in ([55]) we have chosen the ground state of the effective 
Hamiltonian ([M)) but this state has been shifted so that 
the density maximum is not localized around ip ~ but 
around ip — t:. In the present case the delocalization of 
the solitons is stronger than in the case when the solitons 
attract each other. Moreover, initial repulsion of the soli- 
tons weakens in time, i.e. the probability density K(^{p) 
becomes nearly uniformly distributed with small maxima 
around both ip ^ n and p = 0. 

In order to analyze the delocalization of the solitons 
quantitatively we present time evolution of the inverse 
participation ratio (IPR) in Fig. [5] The inverse partici- 
pation ratio. 



IPR = 




is the length scale on which the center of mass of a soliton 
is localized [i^ . While the delocalization of the solitons 
in the case of the attractive soliton interactions is not 
very big (Fig. [S^-b) , the solitons that initially repel each 
other reveal substantial delocalization on the time scale 
presented in Fig. [S] Indeed, the IPR shown in Fig. |BJ:-d 
is comparable to the size of the system L = &. 

In Ref. [i^, the many-body evolution of a BEC, where 
all atoms occupy a double soliton solution of the GPE, 
has been analyzed by numerical simulations. The authors 
consider the cases where the relative phase between the 
solitons is precisely defined and it is either or tt. In the 
both cases they observe a loss of phase coherence between 
the solitons which leads to dramatic fragmentation of the 
condensate. Our results show that if the solitons attract 
each other and the initial state for the phase degree of 
freedom is chosen properly, the phase coherence changes 
very little in the timescale needed for delocalization of 
the soliton centers of mass. Indeed, if we choose as an 
initial state for the phase degree of freedom the ground 
state of the effective Hamiltonian the probability 

density ((36)) remains localized at the bottom of the effec- 
tive potential L^cff, Eq. ([35]). during time evolution. 

If a pair of bright solitons is prepared in a laboratory, 
each density measurement is expected to reveal two soli- 
tons [12, 11^ located at positions with probabilities 
Ki{ri). Averaging atomic density over many realizations 
of the same experiment we obtain a density profile which 
can be compared to the single particle density p{z), see 
Eq. ([^5)) . In Fig. [7] we show the time evolution of p{z). 
Even on a linear scale the delocalization effects can be 
clearly visible, especially in the case when the solitons 
initially repel each other. 

Experimental realization of the Anderson localization 
of a single soliton and observation of the delocalization 
effects in the case of a pair of interacting solitons requires 
the absence of decoherence effects. The most dangerous 
are atomic losses which localize centers of mass of the 
solitons and prevent interference phenomena needed in 
the Anderson localization process. In the present paper 
we concentrate on ^^Rb atoms in the vicinity of the Fesh- 
bach resonance at magnetic field of 155 G. The scattering 
length assumed in our analysis corresponds to the mag- 
netic field of about 169 G. For this magnetic field the 
two- and three-body loss rates measured experimentally 
are lO^^'' cm^ /s and 10~^^ cm^ / s, respectively, and the 
background loss time is 450 s [50|- Assuming that the 
density profile of the atomic cloud is a product of the ID 
soliton density and the density of the harmonic oscillator 
ground state in the transverse directions, the resulting 
lifetime is 160 s, i.e. 5.5 in the units ([2]). This seems 
to be sufficient time to see the delocalization effects an- 
alyzed in the present publication. 
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V. CONCLUSION 

We have considered quantum bright sohtons in the 
presence of an external disorder potentiaL When a single 
soliton is placed in the disorder, Anderson localization of 
its center of mass is predicted by [13] ■ If two solitons 
are present the mutual interaction between them can be 
responsible for the breakdown of the localization. We 
have analyzed this phenomenon within the perturbation 
approach and by means of numerical integration of the 
Schrodinger equation. The perturbation approach shows 
that exponential tails of the probability densities for the 
centers of mass of the solitons disappear due to coupling 
of an unperturbed cigenstate to high energy delocalized 
states that is induced by the solitons interaction. The 
interaction induced delocalization is different from the 
typical situation [8l4Tl| because apart from the relative 



distance between the solitons, the interaction potential 
depends also on the relative phase between them. In the 
numerical integration, where we can afford stronger in- 
teractions than in the perturbation analysis, we observe 
that the localization of the solitons can be totally broken. 

The results obtained in the present publication can be 
verified in experiments. They also indicate that in order 
to observe experimentally the Anderson localization of 
solitons, a single soliton has to be excited and preparation 
of soliton trains should be avoided. 
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